METHOD AND APPARATUS FOR PROVIDING MUTUAL TIME 



DIFFERENCE DETERMINATION OF BASE STATION SIGNALS IN A 
CELLULAR COMMUNICATION SYSTEM 

5 PRIORITY 

This application claims priority to an application entitled "Method of 
mutual time difference determination of base station signals in cellular 
communication system" filed in the Russian Industrial Property Office on 
November 15, 2002 and assigned Serial No. RU2002 130594, the contents of 
10 which are incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

15 The present invention relates to a radio communication system, in 

particular to a method and an apparatus for providing mutual time difference 
determination of base station signals in a cellular communication system. 

2. Description of the Related Art 

Third Generation (3G) CDMA cellular base stations are non-synchronous 
20 in for example the Third Generation Partnership Project (3GPP) Frequency 
Division Duplex (FDD) mode. Here, the synchronism of base stations (BS) 
refers to the synchronization of the downlink transmission signals of different 
base stations. 

Mutual time difference determination is required to determine the mobile 
25 users' locations and to reduce time and hardware resources at the initial base 
station signal search by mobile stations and to decrease the stored data amount in 
the course of the soft handoff procedure. The mutual time difference 
determination during the soft handoff procedure is described in the Russian 
patent #2137314, Int. cl.6 H 04 L 27/30, which is incorporated herein by 
30 reference. 



One of the prior art solutions described in "Synergies Between Satellite 
Navigation and Location Services of Terrestrial Mobile Communication", G. 
Hein, B. Eissfeller, V. Oehler, Jon O. Winkel, Institute of Geodesy and 
Navigation, University FAF Munich, ION GPS 2000, 19-22 September 2000, 
5 Salt Lake City, UT, which is incorporated herein by reference proposes the 
application of location measurement units receiving BS signals and defining the 
BS their mutual time difference to determine mutual time differences of base 
stations. 

However, the disclosed solution described above has the following 
10 problems: 

First, non-line of sight multipath signal propagation from base stations to 
location measurement units causes low accuracy when determining mutual time 
differences of base station signals. 

Second, it is impossible to get the mutual time difference of the base 
15 station signals when the direct measurement is unobtainable. 

The prior art solution to overcome such problems is the solution 
described in 3 GPP Technical Specification (TS) 25.305 V3.7.0, 2001 12, "Stage 
2 Functional Specification of UE Positioning in UTRAN", which is incorporated 
herein by reference. 

20 This solution proposes to receive base station signals, from which the 

mutual time difference is determined by the location measurement unit which is 
located in a position with known coordinates. Fig. 1 helps in understanding the 
method of mutual time difference determination of base station signals in a 
cellular communication system performed by 3GPP TS 25.305 V3.7.0, 2001 12. 

25 Fig. 1 shows base stations 1 and 2, the mutual time difference of which is 

to be determined, a location measurement unit (LMU) 3, base station controller 
(BSC) 4 and a mobile user location center (MULC) 5. 

The LMU 3 receives a signal of the base station (BS) 1 and a signal of 
the BS 2 and executes a specified number of sequential measurements of the 

30 signals' time difference. Further, the LMU 3 averages the sequential 
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measurements of time difference of BS 1 and BS 2 signals, thus obtaining the 
averaged measured time difference of the signals of the given base stations. 

The LMU 3 determines the accuracy of the averaged measured time 
difference, for example, by the signal-to-noise ratio of the BS 1 and BS 2. In this 
5 case a value linearly connected with an error of the averaged measured time 
difference to the true value of the time difference of the BS 1 and 2 signals 
received by the location measurement unit 3 is selected as the averaged measured 
time difference accuracy. 

The averaged measured time difference of the BS 1 and BS 2 signals and 
10 its accuracy are transmitted from the LMU 3 to the BS 1 or BS 2 using the 
current radio interface and then transmitted from the respective BS to the base 
station controller 4 by a wire communication line connecting the BSC 4 and the 
base station. 

The BSC 4 determines the mutual time difference of the BS 1 and BS 2 
15 signals by the averaged measured time difference of the BS 1 and BS 2 signals 
by considering the known mutual location of the BS 1 and BS 2 and the LMU 3. 

The obtained mutual time difference of the BS 1 and BS 2 signals and its 
accuracy are transmitted from the BSC 4 to the mobile user location center 5 for 
further application. 

20 Thus, the following main features of implementing the method can be 

emphasized from the description of the prior art method of mutual time 
difference determination of base station signals in a cellular communication 
system: 

In each location, a measurement unit performs the sequential 
25 measurement of mutual time difference of signals of at least two base stations, 
signals of which are received in the location measurement unit, averaging of 
these time difference measurements results in the averaged measured mutual time 
difference of these base stations, and determines its accuracy. 

The averaged measured time differences and their accuracies are 
30 transmitted from each location measurement unit to one of the base stations 
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which signals are received by the location measurement unit, and then to the base 
station controller that controls the base station. 

The mutual time difference of signals is determined for each base station 
pair for the base station controller using the averaged measured time difference 
5 of the given base station signals. 

However, the known method has a number of important disadvantages. 

First, determination of the mutual time difference of the base station 
signals by the averaged measured time difference of these base station signals 
may be insignificantly accurate. The inaccuracy is caused by the fact that the 
10 estimate of the mutual time difference of the base station signals is subject to the 
impact of noise errors, intra-system interference and multipath errors. 

The difference between the delays of line of sight base station signals 
propagation to the location measurement unit can be removed with the use of the 
known coordinates of the base stations and location measurement unit. 

15 When A ^-> 2 is marked as an estimate of the mutual time difference of 

signals of the base station pair, especially the signals of the first BS to the signal 

of the second BS, and A ^->2 is marked as the true value of the mutual time 
difference of signals of the base stations, the difference between the estimate 

^- >2 and the true value A ^->2 can be expressed by the following Equation 1 . 
20 Equation 1 

^l->2 ~ ^]_>2 = € noise + £ multipath, \ ~ £ multipath ,2 

In equation 1, noise is an error defined by noise and intra-system 

interference, multi P ath > i is a multipath error of the first BS, equal to the difference 
between the real propagation time of the first BS signal to the location 
25 measurement unit and the known propagation time of the first BS signal to the 

location measurement unit, and m,,// ^ rt ' 2 is a multipath error of the second BS, 
equal to the difference between the real propagation time of the second BS signal 
to the location measurement unit and the known propagation time of the second 
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BS signal to the location measurement unit. 

The averaging of the estimate A ^->2 of the mutual time difference of 
signals of the first BS and second BS in the location measurement unit results in 

the reduction of the noise error value noise . The multipath error difference 

5 s muiti P ath,\ ~ e muMpa*,?. j s invariable since it is defined by the mutual location of the 
first and second BS and the location measurement unit and the surrounding 
scattering objects (buildings, mountains, hills, etc.) as well. 

Hence, the accuracy of the claimed method of mutual time difference 
determination of base station signals in a cellular communication system may be 
10 insufficient for a location. 

Second, there may be a situation when there is no direct time signal 
difference measurement between any base stations, even when their mutual time 
difference is required. 

Fig. 2 illustrates the above situation, demonstrating the base stations 6, 7 
15 and 8, location measurement units 9 and 10 and a building 11. Each of the BS 6, 
BS 7 and BS 8 transmits the first and second signals which comprise a group 
signal. 

The base station group signal refers to a signal transmitted from the base 
station and having the Synchronization CHannel (SCH), Common Pilot Channel 
20 (CPICH) and Primary Common Control Physical CHannel (P-CCPCH) and other 
channels as well. 

Each of the LMU 9 and the LMU 10 receives signals transmitted from 
surrounding base stations and measures the signals time difference. 

Referring to FIG. 2, the LMU 9 receives the first signals of the BS 6 and 
25 BS 7 and measures their time difference. The LMU 10 receives the second 
signal of the BS 7 and the first signal of the BS 8 and measures their time 
difference. 

However, the second signal of the BS 6 is held by the building 11 and 
cannot be received by the LMU 10. The second signal of the BS 8 is also held by 
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the building 1 1 and cannot be received by the LMU 9. 

When the direct measurement of the time difference of the signals of the 
BS 6 and BS 8 cannot be performed as described above, the conventional 
methods do not propose any solution capable of determining the mutual time 
5 difference between the signals of the BS 6 and BS 8. 

SUMMARY OF THE INVENTION 

Accordingly, the present invention has been made to solve the above- 
mentioned problems occurring in the prior art, and an object of the present 

10 invention is to provide a method and an apparatus which can determine a mutual 
time difference of signals transmitted from base stations even when the signals 
are blocked by a certain object. 

It is another object of the present invention is to provide a method and an 
apparatus which can improve the determination accuracy of the mutual time 

15 difference of base station signals. 

In order to substantially accomplish this object, there is provided a 
method for determining mutual time difference between base station signals in an 
asynchronous code division multiplexing access (CDMA) system, the method 
comprising the steps of: measuring mutual time difference of signals transmitted 

20 between at least two base stations; determining all possible paths between said at 
least two base stations; and providing weights to the measured mutual time 
difference for said all possible paths. 

In accordance with another aspect of the present invention, there is 
provided an apparatus for determining mutual time difference between base 

25 station signals in an asynchronous code division multiplexing access (CDMA) 
system, the apparatus comprising: a location measurement unit for measuring 
mutual time difference of signals transmitted between at least two base stations; a 
mobile user location center for receiving the mutual time difference of the signals 
measured by the location measurement unit, determining all possible paths 

30 between said at least two base stations, and providing weights to the measured 
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mutual time difference for said all possible paths. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 The above and other objects, features and advantages of the present 

invention will be more apparent from the following detailed description taken in 
conjunction with the accompanying drawings, in which: 

Fig. 1 is a block diagram illustrating an example of a positioning system 
according to a Third Generation Partnership Project Technical Specification 
10 (3GPP TS) 25.305 V3.7.0, 2001 12, "Stage 2 Functional Specification of User 
Equipment (UE) Positioning in Universal Mobile Telecommunications System 
Terrestrial Radio Access Network (UTRAN)"; 

Fig. 2 is a block diagram illustrating an example of an arrangement of 
base stations, location measurement units and buildings when there is no direct 
15 measurement of the time difference of base station signals; 

Fig. 3 is a block diagram illustrating an example of a radio 
communication system in accordance with an embodiment of the present 
invention; 

Fig. 4 is a block diagram illustrating an example of a location 
20 measurement unit for determining the averaged measured time difference of base 
station signals in accordance with an embodiment of the present invention; 

Fig. 5 is a block diagram illustrating an example of a radio 
communication system having base stations as the vertexes and the adjusted time 
differences of the base station signals as the arcs in accordance with an 
25 embodiment of the present invention; 

Fig. 6 is a block diagram illustrating an example of the forming of paths 
from the first base station to the second base station for determining mutual time 
difference of signals in accordance with an embodiment of the present invention; 

Fig. 7 is a flow diagram illustrating an example of a base station in 
30 accordance with an embodiment of the present invention; 
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Fig. 8 is a flow diagram illustrating an example of the base station 
operation in accordance with an embodiment of the present invention; 

Fig. 9 is a flow diagram illustrating an example of a method of operation 
for the base station controller in accordance with an embodiment of the present 
5 invention; 

Fig. 10 is a flow diagram illustrating another example of a method of 
operation for the base station controller in accordance with an embodiment of the 
present invention; 

Fig. 11 is a flow diagram illustrating an example of a method of 
10 operation for the mobile user location center in accordance with an embodiment 
of the present invention; 

Fig. 12 is a flow diagram illustrating an example of a method for 
performing a mobile user location procedure in accordance with an embodiment 
of the present invention; and 
15 Fig. 13 is a flow diagram illustrating an example of a method for 

performing accurate base station location determination in accordance with an 
embodiment of the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

20 

Embodiments of the present invention will be described with reference to 
the accompanying drawings. In the following description of the present 
invention, a detailed description of known functions and configurations 
incorporated herein will be omitted for conciseness. 

25 In acquiring more precise time difference information between base 

stations for measurement of locations in an asynchronous mobile communication 
system, in which no synchronization is made between the mobile stations, the 
embodiment of the present invention provides a method that is capable of 
utilizing mutual time difference information between various base stations, so as 

30 to reduce error in the time difference information which may be caused when a 
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first and second base stations forming a pair among the various base stations are 
located in a downtown area in which propagation delay may frequently happen. 
It should be noted that the embodiment of the present is not limited to base 
stations in downtown areas, the base stations can be used in an environment 
5 where signals can be blocked or affected by physical objects. 

In order to help in understanding the embodiment of the present 
invention, an example of the method of mutual time difference determination of 
base station signals in a cellular communication system proposed by the present 
invention will be described with reference to Figs. 3 to 6. 
10 Fig. 3 is a block diagram illustrating an example of a radio 

communication system in accordance with an embodiment of the present 
invention. 

Referring to FIG. 3, the cellular radio communication system includes a 
plurality of base stations (BSs) 12, 13, 14, 15 and 16, location measurement units 

15 (LMUs) 17, 18, 19 and 20 for receiving at least two base station signals from the 
BSs 12, 13, 14, 15 and 16 through a radio interface, a base station controller 
(BSC) 21 connected through wire communication lines to the BSs 12, 13, 14, 15 
and 16 to control the BSs 12, 13, 14, 15 and 16, and a mobile user location center 
(MULC) 22 connected to the BSC 21 through a wire communication line. The 

20 (MULQ22 could be embedded in the BSC 21 . 

The BS 12 transmits a signal being its group signal. The BS 13 transmits 
the first and second signals being its group signal. The BS 14 transmits the first 
and second signals being its group signal. The BS 15 transmits the first and 
second signals being its group signal. The BS 16 transmits the first and second 

25 signals being its group signal. 

The base station group signal refers to a signal transmitted from the base 
station comprising the following channels: Synchronization CHannel (SCH), 
Common Pilot Channel (CPICH), Primary Common Control Physical CHannel 
(P-CCPCH) and other channels as well. 

30 The LMU 17 receives the signal of the BS 12 and the first signal of the 
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BS 13, and executes sequential measurements of the time difference of the first 
signal of the BS 13 relating to the signal of the BS 12. These time difference 
measurements are averaged, thus obtaining the averaged measured time 

difference A ' 13 -> 12 - 17 of the first signal of the BS 13 regarding the BS 12 signal. Its 

5 accuracy °" 13 -> 12 ' 17 is determined. 

The LMU 17 also receives the first signal of the BS 13 and the first 
signal of the BS 14, and executes sequential measurements of the time difference 
of the first signal of the BS 13 relating to the signal of the BS 14. These time 
difference measurements are averaged, thus obtaining the averaged measured 

10 time difference A ' 13 -> 14 > 17 of the first signal of the BS 13 with regard to the first 

signal of the BS 14. Its accuracy °" 13 -> 14 > 17 is determined. 

The LMU 17 also receives the first signal of the BS 12 and the first 
signal of the BS 14, and performs sequential measurements of the time difference 
of the BS 12 signal with respect to the first signal of the BS 14. These time 
15 difference measurements are averaged, thus obtaining the averaged measured 

time difference A * 12 -> 14 - 17 of the BS 12 signal with respect to the first signal of the 

BS 14. Its accuracy a u->w is determined. 

The LMU 20 receives the first signals of the BS 15 and BS 16, and 
performs sequential measurements of the time difference of the first signal of the 
20 BS 16 with respect to the first signal of the BS 15. These time difference 
measurements are averaged, thus obtaining the averaged measured time 

difference A ' 16 ->i5> 2 ° of the first signal of the BS 16 with respect to the first signal 

of the BS 15. Its accuracy cri6 -> 1520 is determined. 

The LMU 18 receives the second signals of the BS 13 and BS 15, 
25 executes sequential measurements of the time difference of the second signal of 
the BS 15 with respect to the second signal of the BS 13. These time difference 
measurements are averaged, thus obtaining the averaged measured time 
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difference >s->im« of the second signal of the BS 15 with respect to the second 

signal of the BS 13 5 and its accuracy cri5 -> 1318 is determined. 

The LMU 19 receives the second signals of the BS 14 and BS 16, 
performs sequential measurements of the time difference of the second signal of 
5 the BS 14 with respect to the second signal of the BS 16, averages these time 
difference measurements, thus obtaining the averaged measured time difference 

A'i4->i6,i9 Q f t ^ e secon( j signal of the BS 14 with respect to the second signal of 

the BS 16. Its accuracy CTl4 -> 1619 is determined. 

To illustrate, a process of obtaining the averaged time difference 

10 A ' 14 - >16 < 19 of the second signal of the BS 14 with respect to the second signal of 

the BS 16 and determining its accuracy °i4->i6,i9 w jij ^ e explained below with 
reference to Fig. 4. 

Herein, the radio communication system is assumed to be the 3 GPP 
system in the Frequency Division Duplex (FDD) mode. Base stations of the 

15 system transmit signals containing SCH, CPICH and P-CCPCH. 

Fig. 4 is a block diagram illustrating an example of a location 
measurement unit for determining the averaged measured time difference of base 
station signals in accordance with an embodiment of the present invention. 

Referring to Fig. 4, the location measurement unit for determination of 

20 the averaged time difference of signals of the BS 14 and BS 16 comprises an 
antenna 23, analog receiver 24, searcher 25 of the BS 14 group signal 
(hereinafter referred to as searcher 25), decoder 26 of the Primary Common 
Control Physical Channel of the BS 14 (hereinafter referred to as decoder 26), 
searcher 27 of the BS 16 group signal (hereinafter referred to as searcher 27), 

25 decoder 28 of the Primary Common Control Physical Channel of the BS 16 
(hereinafter referred to as decoder 28) and a unit 29 of determination of the 
averaged measured time difference of signals of the BS 14 and 16 (hereinafter 
referred to as unit 29), wherein the input of the antenna 23 is the input of the 
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location measurement unit, the output of the antenna 23 is connected to the input 
of the analog receiver 24, the output of which is linked to the inputs of the 
searcher 25 and the searcher 27 and to the first inputs of the decoder 26 and 
decoder 28, the output of the searcher 25 is connected to the first input of the unit 
5 29 and to the second input of the decoder 26, the output of which is connected to 
the second input of the unit 29, the output of the searcher 27 is connected to the 
third input of the unit 29 and to the second input of the decoder 28, the output of 
which is connected to the fourth input of the unit 29, the output of which is the 
output of the location measurement unit. 

10 In this case, the analog receiver can be embodied in the manner provided 

in the US Patent US No. 5,103,459 titled "System and Method for Generating 
Signal Waveforms in a CDMA Cellular Telephone System", Int. C1.5 H 04 L 
27/30, the contents of which is incorporated herein by reference. 

The searcher 25 and the searcher 27 can be embodied in the manner 

15 described in 3 GPP TS 25.214 V3.9.0 (2001-12), Physical layer procedures 
(FDD), Annex C: Cell search procedure and in "Cell Search in W-CDMA", Yi- 
Pin Eric Wang and Tony Ottosson, IEEE JOURNAL ON SELECTED AREAS 
IN COMMUNICATIONS, VOL. 18, NO. 8, AUGUST 2000, the contents of 
which is incorporated herein by reference. 

20 The decoder 26 and the decoder 28 can be embodied like the coherent 

RAKE receivers of Sadayuki ABETA, Mamoru SAWAHASHI, and Fumiyuki 
ADACHI, "Performance Comparison between Time-Miltiplexed Pilot Channel 
and Parallel Pilot Channel for Coherent Rake Combining in DS-CDMA Mobile 
Radio", IEICE Trans. Commun., Vol. E81-B, No. 7, July 1998, the contents of 

25 which is incorporated herein by reference. 

An input signal of the location measurement unit containing BS 14 and 
BS 16 group signals is applied to the input of the antenna 23 and is applied from 
the output of antenna 23 to the input of the analog receiver 24. The signal is 
transmitted from the output of the analog receiver 24 to the inputs of the searcher 

30 25 and searcher 27 and to the first inputs of the decoder 26 and the decoder 28. 
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The searcher 25 searches the group signal of the BS 14 in the searching window 
using 38400 chips with the Synchronization Channel and Common Pilot 

Channel. A signal of the BS 14 is assumed to be acquired in the position P K 
Additionally, the searcher 25 determines the number of the primary scrambling 
5 code of the BS 14. A signal containing the values of the acquired time position 

^ of the BS 14 signal is applied from the output of the searcher 25 to the first 
input of the unit 29 and to the second input of the decoder 26. A signal containing 
the number of the primary scrambling code of the BS 14 is applied from the 
output of the searcher 25 to the second input of the decoder 26. With the 

10 obtained value of the acquired time position r i of the BS 14 signal and the 
obtained number of the primary scrambling code of the BS 14 the decoder 26 de- 
scrambles, demodulates and decodes the Primary Common Control Physical 
Channel, thus obtaining the value of the System Frame Number (SFN) of the BS 
14 at the moment of the transmission of the first chip of this frame. This SFN 

15 value of the BS 14 is labeled by ^ FN \ . A signal containing the determined value 

SFN * of the SFN of the BS 14 is applied from the output of the decoder 26 to the 
second input of the unit 29. 

The searcher 27 searches the group signal of the BS 16 within the 
searching window using 38400 chips with the Synchronization Channel and 

20 Common Pilot Channel. It is assumed that the BS 16 signal is acquired in the 

p 

position 2 . Additionally, the searcher 27 determines the number of the primary 

scrambling code of the BS 16. A signal containing the value of the acquired time 
p 

position 2 of the BS 16 signal is applied from the output of the searcher 27 to 

the third input of the unit 29 and to the second input of the decoder 28. A signal 

25 containing the number of the primary scrambling code of the BS 16 is applied 

from the output of the searcher 27 to the second input of the decoder 28. With 

p 

the obtained value of the acquired time position 2 of the BS 16 signal and the 
obtained number of the primary scrambling code of the BS 16 the decoder 28 de- 
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scrambles, demodulates and decodes the Primary Common Control Physical 
Channel, thus obtaining the value of the SFN of the BS 16 at the moment of the 
transmission of the first chip of this frame. This SFN value of the BS 16 is 

labeled by SFN 2 m a signal containing the determined value SFN 2 G f the SFN of 
5 the BS 16 is applied from the output of the decoder 28 to the fourth input of the 
unit 29. 

The unit 29 defines the measured time difference At "->w of the BS 14 
signal with respect to the BS 16 signal using the following equation 2. 
Equation 2 

! o A 'i4->i6,i9 = (SFN t - SFN 2 )T fr + {P x -P 2 ) T ch 

T 

In equation 2, fr - length of a frame of 3GPP base station signal of 10 

ms, and Tch - length of one chip of V( 3 10 ) s or around 260 ns. 

The unit 29 averages, for example, some sequentially measured time 
differences of the BS 14 signal with respect to the BS 16 signal obtained by 

15 equation 2, thus obtaining the averaged measured time difference A ' 14 -> 16 > 19 of the 
signal of the BS 14 with respect to the signal of the BS 16. The unit 29 can be 
implemented on the Digital Signal Processor (DSP) using the above algorithm. 

If all location measurement units determine the averaged measured time 
differences of base station signals in a similar manner with the same number of 

20 averages, the accuracies of all averaged measured time differences of base station 
signals can be defined to be the same and equal to 100 ns, for example. 

In a general case with more complex and more accurate determination 
methods of the averaged measured time differences of base station signals or 
with different number of averages in different location measurement units, the 

25 accuracies of different averaged measured time differences of base station signals 
are different and depend on, for example, signal-to-noise ratios of signals, which 
averaged measured time difference to be determined, and/or number of averages. 
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The averaged measured time differences ^^-^w y A* 13 _ >14fl7 ^ Ar 12 _ >1417 

and their accuracies °" 13 -> 12 . 17 ? a n->\4j7 an( j a n~>\4M are transmitted from the LMU 
17 to the BS 13 and then to the BSC 21. 

The averaged measured time difference ^*«*->i3,i8 and its accuracy 

5 °"'5->i3,i8 are transmitted from the LMU 18 to the BS 13 and then to the BSC 21. 

The averaged measured time difference A ' 14 -> 16 > 19 and its accuracy 

^ are transmitted from the LMU 19 to the BS 14 and then to the BSC 21. 

The averaged measured time difference A ' 16 -> 15 > 20 and its accuracy 

0"i6->i5,2o are transmitted from the LMU 20 to the BS 15 and then to the BSC 21. 
10 As will be shown from the following equation 3, the known value of the 

difference of delays at line of sight signal propagation from the BS 13 and the BS 
12 to the location measurement unit 17 is subtracted from the averaged measured 

time difference ^ 13 -> I2 > 17 in the BS controller 21, thus obtaining the adjusted 

time difference ^ 13 - >12 > 17 Q f the signal of the BS 13 with respect to the signal of 
15 theBS 12. 

Equation 3 

^13->12,17 = A*13->12,17 (^~13->17 "" r i2->l?) 

In equation 3, *3->n j s t he signal propagation delay from the BS 13 to 

the location measurement unit 17 and 12_>17 is the signal propagation delay 
20 from the BS 12 to the location measurement unit 17. The signal propagation 
delays can be calculated using the following equation 4. 
Equation 4 

_ V(*i2 -*n) 2 +Qyi2 -y\i) 2 +( z i2 - z nY 
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_ J(xn- X n) 2 +(yn-yn) 2 +{ z n- z nf 

^1 -J_-s.» -7 



M3->17 

c 



In equation 4, c is the light speed, the BS 12 coordinates with 

*i2>.Vi2> z i2 3 t he BS 13 coordinates with ^'^n^ia ? and the location 

„ X V z 

measurement unit 17 coordinates with I7'^i7' 17 -j^e base stations and 

5 location measurement unit coordinates can be determined with the use of the 

Global Positioning System (GPS)and/or Global Navigation Satellite System 

(GNSS ) receiver. 

In this case, the accuracy of the adjusted time difference ^ 13 ~ >12 > 17 i s 
equal to the accuracy of the adjusted averaged measured time difference 

The adjusted time difference A/ »3->i4,i7 of the Q f the BS 13 with 

respect to the signal of the BS 14, the adjusted time difference ^ ,2 - >14 > 17 of the 
signal of the BS 12 with respect to the signal of the BS 14, the adjusted time 

difference ^ 15 - >,3 > 18 of the signal of the BS 15 with respect to the signal of the 
15 BS 13, the adjusted time difference At "->w> of the signal of the BS 14 with 

respect to the signal of the BS 16 and the adjusted time difference ^^ 6 - >15 2 ° of 
the signal of the BS 16 with respect to the signal of the BS 15 are determined in 
the same manner. 

The adjusted time differences A '«3->i2,i7 } M 3 _ >14il7 ^ Af 12 _ >1417 ^ 
20 A A 5 ->i3,. 8> A/ 14 _ >16)19 and Ar i6 _ >l5>20 andtheiraccuracies cr n->i2,n > cr i3->i4,i7 ) 
^i2->i4,i7 ^ C7 15 _ >13 18 ^ cr 14 _ >l6 19 and o- 16 _ >1520 ^g^ttgj from the BSC 2i to 
the MULC 22. 
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Assume, the mutual time difference ^3->u G f t he s ign a i c f the BS 13 
with respect to the signal of the BS 14 should be determined by the MULC 22. 

The true mutual time difference value of the signal of the BS 13 with 

respect to the signal of the base station 14 is marked by ^3-> 14 9 the true mutual 
5 time difference value of the signal of the BS 13 with respect to the signal of the 

BS 12 is marked by ^ 13 ~ >12 , the true mutual time difference value of the signal 

of the BS 12 with respect to the signal of the BS 14 is marked by A ^2->i4 ? t he 
true mutual time difference value of the signal of the BS 15 with respect to the 

signal of the BS 13 is marked by ^5->i3 ^ the true mutual time difference value 
10 of the signal of the BS 14 with respect to the signal of the BS 16 is marked by 

A ^i4->i6 anc j t he true mutual time difference value of the signal of the BS 16 
with respect to the signal of the BS 15 is marked by ^6->i5 

Since A/ i3->i4 = A? i3->i2 + A 'i2->i4 and 

A 'i3->i4 = - A 'is->i3 " A 'i6->i5 - A 'i4->i6 9 the mutual time difference A *i3->i4 can 
15 be estimated by three methods: 

A *13->14,17 

A ^13->12,17 + A ^12->14,17 H 

~ A *15->13,18 ~~ A ^16->15,20 ~~ A ^4->16,19 

The configuration composed of the BS 12 through BS 16 and the 
20 adjusted time differences A ^->12,17 ? A/ i3->i4,i7 ^ A 'i2->i4,i7 ? A/ 15 _ >13 18 ^ 

A ^i4->i6 s i9 an( j A *i6->15,20 j s p resen ted in the form of a graph depicted in Fig. 5. 

Referring to Fig. 5, the vertexes of the graph are the BS 12 through BS 
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16 and its arcs are the adjusted time differences ^ 13 - >12 > 17 , >14 » 17 9 

A *12->14,17 A ^5->13,18 ^14->16,19 an( j ^16->15,20 

The graph arc directions complying with the directions of the adjusted 
time difference are set. For example, the arc between the BS 13 and BS 14 

5 matches the adjusted time difference A ' 13 -> 14 > 17 of the signal of the BS 13 with 
respect to the signal of the BS 14 and directs from the BS 13 to the BS 14, 
respectively. 

Tree estimates A/ n->i4,i7 ? ^ 13 _ >12>17 + Af 12 _ >14>17 ^ 

-^ 15 _ >13 ,i8-^i6->i5,2o- Af i4->i6,i9 of the mutual time difference A *i3->14 
10 correspond to three paths of the graph from the BS 13 to the BS 14 depicted in 
Fig. 6. 

The three paths depicted in Fig. 6 are the set of all possible paths from 
theBS 13 to the BS 14. 

The vertex of each path is a part of base stations, which are the graph 

15 vertexes. The first vertex of each path is the BS 13 and the last vertex is the BS 
14. The second vertex of the first path is the BS 12; adjacent to the first vertex 
being the BS 13. The second vertex of the third path is the BS 15; adjacent to the 
first vertex being the BS 13. The third vertex of the third path is the BS 16; 
adjacent to the second vertex of the third path being the BS 15. 

20 The adjacent base stations are base stations for which the adjusted time 

difference of their signals has been obtained. For instance, the BS 15 and the BS 

16 are adjacent, since the adjusted time difference ^ 16 - >15 > 20 G f their signals has 
been obtained. 

The direction of each path is defined from the BS 13 to the BS 14. 
25 Six adjusted time differences have been obtained. Thus, the digits from 

1 to 6 number these adjusted time differences by the following equation 5. 
Equation 5 
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Af, 




>12,17 


At 2 




->14,17 


At 3 


= AF,, 


->14,17 


& 4 


= Af 15 . 


->13,18 


At 5 


= Af 14 . 


->16,19 

5 


A' 6 


= A? 16 


->15,20 



The terms of the vector of the adjusted time differences and the vector of 
the adjusted time difference accuracies are here introduced. 

The vector of the adjusted time differences having the length of six, 
10 being number of the obtained adjusted time differences, means a vector 
expressed by the following equation 6. 

Equation 6 



■a/," 




A^13->12,17 


At 2 




A 'l3->14,17 


At 3 




A ^12->14,17 


At, 




A*15->13,18 


At s 




A ^14->16,19 


_At 6 _ 




_A*16->15,20_ 



In equation 6, ^-th element tp of the vector A ' , where p gets the 
15 values from 1 to 6 , equals to the p -th adjusted time difference of the signals of 
the base stations. 

The vector of the adjusted time difference accuracies ° having the 
length of six, being number of the obtained adjusted time differences, refers to a 
vector expressed by the following equation 7. 
20 Equation 7 
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a = 













°"l3->14,17 






°12->14,17 






^"lS->13,18 






Cr i4->16,I9 






_ Cr i6->15,20_ 



In equation 7, ^ -th element of the vector ^ , where ^ gets the values 
from 1 to 6 , equals to the accuracy of the p -th adjusted time difference of the 
base station signals. 

5 The path vector enumerating the adjusted time differences of the base 

station signals comprising this path is formed for each path from the formed set 
of all possible paths from the BS 13 to the BS 14. 

The length of each path vector of six, being number of the adjusted time 
differences, is set. 

10 The ^-th element of the w -the path vector (where p takes the values 

within 1 - 6 and u takes the values from 1 to 3 ) is equal to 1 if the path has an 
arc corresponding to the p -th adjusted time difference and the direction of the 
pass of the arc of the path complies with the direction of the p -th adjusted time 

difference, where the direction of the p -th adjusted time difference of the lp -th 

1 5 base station with respect to the Jp -th base station is defined from the lp -th base 

station to the ^-th base station; it is equal to -1 if the path has an arc 
corresponding to the p -th adjusted time difference and the direction of the pass 
of the arc of the path is reverse to the direction of the p -th adjusted time 
difference; otherwise it is equal 0 . 
20 The above described action results in three path vectors B l ? B 2 and B 3 

corresponding to three paths from the BS 13 to the BS 14. These path vectors 
equal to the following equation 8. 
Equation 8 
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0 




-1 


0 




0 




-1 


0 




0 




-1 



The metric of the u -th path vector where u takes the values from 1 to 3 
is determined by the following equation 9. 
Equation 9 

5 >=' 

In equation 9, p is p -th element of the u -th path vector, and p is p - 
th element of the vector <* of the adjusted time difference accuracies of the base 
stations. 

— 2 2 

Thus, the metric of the first path vector B x i s a \ + °" 3 , the metric of the 

- 2 2 ~ 3 

10 second path vector ^ is 2 ; and the metric of the third path vector B is 

A group of path vectors is selected for each pair of base stations from the 
set of all possible path vectors of this pair of base stations, wherein the selected 
group of the path vectors contains each of the obtained adjusted time differences; 
1 5 wherein the number of applications of each obtained adjusted time difference of 
the selected group of path vectors must not exceed the number of applications of 
this adjusted time difference of any other group of path vectors obtained from the 
set of all possible path vectors, wherein values of path vector metrics of the 
selected group must not exceed values of path vector metrics of any other path 
20 vector group obtained from the set of all possible path vectors. 

These three criteria of the path vector group selection will be explained 
in greater details further. 

The selected path vector group of this example complies with the set of 
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all possible path vectors from the BS 13 to the BS 14. 

Three estimates of the mutual time difference ^*i3->w of the signal of the 
BS 13 with respect to the signal of the BS 14 can be formed using the vector of 
the adjusted time differences A * and three formed path vectors B l 9 B 2 and B 3 . 
5 In this case, the r-th estimate of the mutual time difference A ^3->i4 where r takes 
the values from 1 to 3 is given by the following equation 10. 

Equation 10 

In equation 10, p is the p -th element of the r -th vector of the path # r , 

10 and ^ tp is the ^-th adjusted time difference (^-th element of the mutual time 
difference vector A* ). 

An error r of the ^-th mutual time difference estimate ,3 ~ >14 with 

respect to the true value of the time difference A ^3->i4 Q f the signal of the BS 13 
with respect to the signal of the BS 14 is expressed by the following equation 11. 
15 Equation 11 

3 

p=l 

The correlation matrix between the estimate errors of the mutual time 

difference A ^3->i4 obtained by separate path vectors is formed with the use of the 
vector of accuracies of the adjusted time differences of the signals of the base 
20 station & and three generated path vectors B l 9 B 2 and # 3 , and the correlation 

matrix size is 

[3x3] 

The correlation matrix element having index marks r * and 

r 2 (where ^ and Vl take the values from 1 to 3 ) is expressed by the following 
equation 12. 

Equation 12 
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r, ,r 2 L r \ r 2 J ^— i P P P 

p=l 

In equation 12, L n» J j s the correlation coefficient between the 



error 



At. 



r > of the 1 -th estimate of the mutual time difference 13 - >14 and the error r * of 

r At b* x n 

the 2 -th estimate of the mutual time difference ,3 ~ >14 , p is the p -th element of 

5 the Vx -th path vector B h , p is the ^ -th element of the r * -th path vector B' 2 . 

The correlation matrix 13 ~ >14 of this example is expressed by the 
following equation 13. 
Equation 13 

- 2 0 0 



13->14 



0 
0 



a\ 0 
0 ct 2 + cr 5 2 + a\ 
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In general, the correlation matrix is non-diagonal. 
W 

The matrix r » 3 ->i 4 , being inverse to the obtained correlation matrix 
^i3->i4 i s formed. The matrix W ™->™ inverse to the formed correlation matrix 



B->i4 0 f the present example is expressed by the following equation 14. 
Equation 14 

1 



^13->14 ( ^13->14 ) ~~ 
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CT, + <T 3 



0 
0 



1 



0 
0 

1 



<7 4 2 +o- 5 2 +cx 6 2 



With the use of the formed path vector group consisting of the path 
vectors B 1 t B 2 and B 3 and the formed matrix ^»->m inverse to the generated 

i?- 

correlation matrix l! - >M , six weights of the adjusted time differences are 
produced, and the weight of the p -th adjusted time difference a " (where p 
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takes the values from 1 to 6) may be calculated by the following equation 15. 
Equation 15 



3 3 



(*;+*;) 



P 3 3 

r^l r 2 =l 

In equation 15, is an element of the formed matrix W ^->^ with the 

5 indices 1 and 2 , where 1 and 2 have the values from 1 to 3 , p is the "-th 

element of the -th path vector £ r ' , and p is the p -th element of the ^ 2 -th path 
vector . 

The weights of the adjusted time differences of the present example 
equal to those shown in the following equation 16. 
10 Equation 16 

a\{a\ +o- 5 2 +C7 6 2 ) + C7 2 (cr 2 + 0 - 2 ) + (cr 1 2 + cr 2 )(a 2 +a 2 +cr 2 ) 



The mutual time difference A ^3->i4 Q f t h e signal of the BS 13 with 
15 respect to the signal of the BS 14 is determined by the following equation 17. 
Equation 17 



6 

A 'l3->14 = IX A 'p 



That is, it is determined as the weighted sum of all adjusted time 

differences ^ tp of the signals of the base stations, wherein the weights Up of the 
20 adjusted time differences of the signals of the base stations, where p has values 
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from 1 to 6 9 are used as weights. 

The radio communication cellular system comprising the base station 
controllers, base stations, location measurement units and a mobile user location 
center, which can practice the embodiment of the present invention as described 
5 above, will now be discussed, wherein each base station is controlled by one base 
station controller, each location measurement unit receives signals of at least two 
base stations; and a signal of each base station is received by at least one location 
measurement unit. 

The base stations, base station controllers and the mobile user location 
10 center can be embodied in a manner described in Patent WO #99/57826 titled 
"Method Of Synchronization Of A Base Station Network", 4 May 1998, Int. CI. 
H 04 J 3/06, H 04 B 7/26, the contents of which is incorporated herein by 
reference. 

Below is the brief description of units to include in the listed components 
15 of the cellular radio communication network to embody the method of mutual 
time difference determination of base station signals in a cellular communication 
system in accordance with an embodiment of the present invention. 

Each base station transmits a signal being its group signal. 

Each location measurement unit measures sequentially the time 
20 difference of signals of at least two base stations, which are signals to be received 
by the location measurement unit, these time difference measurements are 
averaged, thus obtaining the averaged measured time difference of signals of the 
base stations. The accuracy of the averaged measured time difference of signals 
of the base stations is determined. The averaged measured time differences and 
25 their accuracies are transmitted from each location measurement unit to one of 
the base stations. The signals are received by the location measurement unit, and 
then to the base station controller controlling the base station. 

Referring to Fig. 7, the BS 31 should comprise at least a receiver 34 of a 
signal containing the averaged measured time difference of signals of a base 
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station pair and its accuracy (hereinafter referred to as receiver 34), a transmitter 
35 of a signal containing the averaged measured time difference of signals of a 
base stations pair and its accuracy (hereinafter referred to as transmitter 35). 
With reference to Fig. 8, the BS 31 executes the following steps. 
5 In step S25, the receiver 34 receives a signal, containing the averaged 

measured time difference of signals of a base station pair and its accuracy, 
transmitted from the LMU 30. In step S30, the transmitter 35 transmits a signal 
containing the averaged measured time difference of signals of a base station pair 
and its accuracy to the BSC 32. 
10 Referring to Fig. 9, the BSC 32 should comprise at least a receiver 36 of 

a signal containing the averaged measured time difference of signals of a base 
stations pair and its accuracy, a computing unit 37 and a transmitter 38 of a signal 
containing the adjusted time difference of signals of a base station pair and its 
accuracy. 

15 As shown in Fig. 10, the following operations are performed using the 

BSC 32. In step 130, the receiver 36 receives a signal containing the averaged 
measured time difference of signals of a base station pair and its accuracy 
transmitted by the BS 3 1 . In step 135, the computing unit 37 subtracts the known 
value of the delay difference at the line of sight signal propagation from the first 

20 BS and second BS of the pair to the location measurement unit where this 
averaged measured time difference has been obtained from the averaged 
measured time difference signals of each pair of base stations, thus obtaining the 
adjusted time difference of signals of this base stations pair. In step 140, the 
transmitter 38 transmits a signal containing the adjusted time difference of the 

25 base stations pair and its accuracy to the MULC 33. 

Referring to Fig. 11, the MULC 33 comprises at least a receiver 39 of a 
signal containing the adjusted time difference of signals of the base stations pair 
and its accuracy (hereinafter referred to as the receiver 39) and a computing unit 
40. 

30 The following steps are performed by the mobile user location center as 
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shown in Fig. 12. In step 240, the receiver 39 receives a signal containing the 
adjusted time difference of signals of the base station pair and its accuracy 
transmitted by the BSC 32. 

In step 245, the computing unit 40 forms a set of all possible paths from 
5 the first BS to the second BS for each pair of base stations. Said all possible paths 
refer to paths between base stations adjacent to a terminal which is an object of 
the measurement, and the paths between the base stations adjacent to the terminal 
include non-line-of-sight multipaths. Assume, the radio communication cellular 
system consists of L base stations. All base stations are assumed to be numbered 
10 from 1 through L . 

Below is the explanation of the term "adjacent base station". The BS i is 
adjacent to the BS J , (wherein 1 and J takes the values from 1 to if at least 
one adjusted time difference of a signal of the BS 1 is obtained with respect to a 
signal of the BS J . If p adjusted time differences of a signal of the BS Ms 
15 obtained with respect to a signal of the BS J , these base stations are adjacent to 
each other p times. 

The set of adjacent base stations is formed for each base station in the 
following manner. Assume, the number of the base stations adjacent to the BS i 
is Qi , where i takes the values from 1 to I. The set of base stations adjacent to 

20 the BS i is defined by the following equation 18. 
Equation 18 

A = {W 4) 

In equation 18, q. r -th element d l of the set D. , where q. takes the values 
from 1 to Q. , is the number of the q i -th base station adjacent to the BS i . If any 
25 base station is an adjacent base station to the BS i p times, its number is 
included to the D i set p times. L sets D. of the adjacent base stations are 
formed. 

The set of all possible paths from the first BS to the second BS of the 
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pair is formed for each pair of base stations in the following manner. 

Assume, all paths from the BS lfn to the BS Jm have to be determined, 

where lfn and Jm have the values from 1 to L . This is performed by the 

sequential selection of the following: all base stations adjacent to the BS lfn , all 

5 base stations adjacent to the base stations which are adjacent to the BS lm , all 
base station adjacent to the base stations of adjacent base stations which are 

adjacent to the BS lfn , etc. 

The sequential selection is performed: the selection sequence to be 

id* d* d* \ 

determined uniquely by sequences of the numbers of base station < °' Q >> 

10 of the formed sets ^ of the adjacent base stations; each path to be updated with 
those numbers of base stations, which are not found previously in the sequential 
selection. In general, the lengths of the resulting paths are different and the 
length of each path is less or equal to number of the base stations L . Each path 
is a sequence of the numbers of the base stations in pass during the sequential 

1 5 selection. The first element of the sequence of base station numbers of each 

formed path is lm . The last element of the sequence of base station numbers of 
each formed path can be either Jm or not. Those formed paths which end at the 
BS Jm are kept, i.e. those formed paths with the last element of the sequence of 
base station numbers of Jm are kept. Number of the paths kept is assumed to be 

20 U m They form a set of all possible paths from the BS ljn to the BS Jm . They are 

numbered with the numbers from 1 to U . 

In step 250, a path vector listing the adjusted time differences of signals 

of base stations included into this path is formed for each path of each formed 

set; the metric of the said path vector is determined. 
25 Assume, P adjusted time differences of the base station signals are 

obtained. They are numbered by the numbers from 1 to P. The terms of an 
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adjusted time difference vector and an adjusted time difference accuracy vector 
are introduced. 

The adjusted time difference vector with the length of P (number of 

the obtained adjusted time differences) is expressed by the following equation 19. 

5 Equation 19 

' At x " 
AL 

At = . 2 
_At p _ 

In equation 19, the p -th element ^ tp of the vector AF (where p takes 
the values from 1 to P ) equals to the p -th adjusted time difference of the base 
station signals. 

10 The vector of the adjusted time difference accuracies ° with the length 

of P (number of the obtained adjusted time differences) is a vector expressed by 
the following equation 20. 
Equation 20 




15 In equation 20, the p -th element (Tp of the vector ° (where p takes the 

values from 1 to P ) equals to the accuracy of the p -th adjusted time difference 
of the base station signals. 

The path vector enumerating the adjusted time differences of the signals 
of the base stations comprised to this path is formed for each path of the 
20 generated set of all possible paths. 

To illustrate, consider this operation for the set of all possible paths from 

the BS l » to the BS J ™ . 

The length of each path vector is set to be P 9 number of the adjusted 
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time differences. 

The P -th element 6f the w -th path vector (where P takes the values 
from 1 to P and u takes the values from 1 to U ) is equal to 1 , if this path has 
an arc corresponding to the p -th adjusted time difference and the direction of the 
5 pass of the arc of the path complies with the direction of the p -th adjusted time 
difference; it is equal to -1 , if the path has an arc corresponding to the p -th 
adjusted time difference and the direction of the pass of the arc of the path is 
reverse to the direction of the P -th adjusted time difference; and otherwise it is 
equal 0. 

10 U vectors of the B x , B u paths from the BS lm to the BS Jm are 

obtained. The formed set of the vectors of all possible paths from the BS l ™ to 

the BS J™ is generally redundant. 

The metric of each formed path vector is determined. For instance, the 

metric of the u -th vector B u of the set of vectors of all possible paths from the 

15 BS l/n to the BS Jm is determined in the manner expressed by the following 
equation 2 1 . 

Equation 21 

p=\ 

In equation 21, p is the p -th element of the u -th vector of the path B u 

20 from the BS l ™ to the BS Jm , and °* is the p -th element of the vector ° of the 

adjusted time difference accuracies. 

In step 255, a group of path vectors is selected for each pair of base 

stations from the set of all possible path vectors of this pair of base stations, 

wherein the selected group of the path vectors contains each of the obtained 
25 adjusted time differences; wherein number of applications of each obtained 

adjusted time difference of the selected group of path vectors must not exceed the 
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number of applications of this adjusted time difference of any other group of path 
vectors obtained from the set of all possible path vectors, wherein values of path 
vector metrics of the selected group must not exceed values of path vector 
metrics of any other path vector group obtained from the set of all possible path 
5 vectors. 

This operation is demonstrated by an example of the selection of a group 
of the path vectors from the set of all possible path vectors from the BS lfn to the 

BS J »> . 

The set of all possible path vectors from the BS lfn to the BS Jm is sorted 
10 in order of increasing path vector metric. The sorting results in the obtaining of 
U vectors of the paths B 1 9 B u 9 where the vector of the path B l has the 
minimum metric and the vector of the path B u has the maximum metric. 

The vector of the adjusted time differences C of the B u path vector is 

defined. The c vector length is P , and the p -th element of the vector c equals 
15 to equation 22 below. 

Equation 22 




That is to say, the element equals to the absolute value of the p -th 
element of the vector of the path B u . 
20 The adjusted time difference vector C of the #\ B Ul ? .... path vectors is 

defined. The length of the vector C is P , and the P -th element of the vector C 
is determined by the manner below: Cp ~~® y if ^ p ^ p otherwise 

To select the path vector group from the set of all possible path vectors 
25 from the BS l » to the BS Jm , B 2 , B u 

are tried to be removed sequentially 
from the path vector set at U -\ steps. In this case the vector of the path B u is 
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attempted to be removed at the u -th step (where u takes the values from 2 to 
^ ) in the following manner: 

the vector C is formed as a vector of the adjusted time differences of the 
vectors of the paths B ul , B u 2 9 B 1 5 
b u = 0 

5 if p for all ^ (where ^ takes the values from 1 to P) for which 

c ^ ~ ® , the vector of the path B u can be removed, 
otherwise, it cannot be removed. 

Assume, we have R path vectors from the BS l,n to the BS Jm . They 
form the selected group of the path vectors B\ B R from the BS to the BS 

10 Jm , which will be further used. 

In step 260, weights of the adjusted time differences of signals of base 
stations are formed for each pair of base stations using the selected group of path 
vectors of these base stations and the obtained accuracies of the adjusted time 
differences of signals of the base stations. That is, for said all paths, weights are 

1 5 provided according to errors of the measured mutual time differences. This 
situation is explained by the forming of the weights of the adjusted time 

difference of BS lfn andBS Jm signals. 

The correlation matrix of errors of mutual time difference of the signal of 

the BS lfn with respect to the signal of the BS Jm , obtained by separate path 
20 vectors is formed using the accuracy vector & of the adjusted time differences of 
the base station signals and the selected path vectors group B\ . . . ? B R m 

The correlation matrix size is [^ x ^]. The correlation matrix element 

with the indices Vx and r * (wherein r * and r * taking the values within 1 - R ) is 
expressed by the following equation 23. 
25 Equation 23 

r \ » r 2 L r l r 2 J -fc-J P P P 
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In equation 23, L^v 6 ^ J j s the correlation coefficient between the -th 
and r2 -th errors of the estimate of the mutual time difference of the signal of the 
BS lfn with respect to the signal of the BS Jm , p is the P -th element of the r * -th 
path vector B r * , and p is the ^ -th element of the r * -th path vector B ri . 

5 The matrix w , inverse to the formed correlation matrix K is produced. 

P weights of the adjusted time differences are generated using the 
selected group of the path vectors consisting of the path vectors B l 9 . .., B R and 
using the generated matrix w , inverse to the formed correlation matrix K , the 

weight of the ^-th adjusted time difference Gp (wherein P taking the values 
10 from 1 to P) being calculated by the following equation 24. 
Equation 24 

n=\ r 2 =l 

CL = — 

u p R R 

In equation 24, Wr ^ is an element of the formed matrix w with the 

indices 1 and 2 , where 1 and 2 take the values from 1 to R , p is the " -th 

15 element of the r * -th path vector B n ? and ^ is the ^ -th element of the r2 -th path 
vector B' 2 . 

In step 265, a mutual time difference of signals of each pair of base 
stations is determined as a weighted sum of all adjusted time differences of 
signals of base stations, wherein the weights of the adjusted time difference of 
20 signals of the base stations formed for the pair of base stations are used as 
weights. 

For example, the mutual time difference ^^J* of the signal of the BS 
%m with respect to the signal of the BS Jm is determined by the following 
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equation 25. 

Equation 25 

p 

'm ~>Jm Z—i P P 

p=] 

That is to say, it is determined as the weighted sum of all adjusted time 
5 differences ^ tp of the base station signals wherein the weights Qp of the adjusted 

time differences of the base station signals formed for the BS lfn and BS Jm , 
where p takes the values within 1 ... P , are used as weights. 

Hereinafter, the joint operation of the cellular system components for the 
implementation of the claimed invention will be described with reference to Fig. 

10 13. Fig. 13 shows the LMU 30, BS 31, BSC 32 and the MULC 33. 

In step 310, the LMU 30 sequentially measures the time difference of 
signals of at least two received base stations. In step 315, the LMU 30 averages 
these time difference measurements, thereby obtaining the averaged measured 
time difference of signals of these base stations. In step 320, the LMU 30 

15 determines the accuracy of the averaged measured time difference of signals of 
these base stations. In step 325, the LMU transmits the first signal containing the 
averaged measured time difference of the signals of the base stations pair and its 
accuracy to the BS 31. In step 330, the BS 31 transmits the second signal 
containing the averaged measured time difference of the signals of the base 

20 stations pair and its accuracy to the BSC 32. 

In step 335, in the BSC 32, the known value of the delay difference at the 
line of sight signal propagation from the first BS and second BS of each pair of 
the base stations to the location measurement unit, where this averaged time 
difference is obtained, is subtracted from the averaged measured time difference 

25 of this pair of the base stations, thus obtaining the adjusted time difference of the 
signals of this base stations pair. 

In step 340, the BSC 32 transmits the third signal containing the adjusted 
time difference of each base stations pair and its accuracy to the MULC 33. 
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In step 345, the MULC 33 forms a set of all possible paths from the first 
BS to the second BS of the pair for each base stations pair. 

Said all possible paths refer to paths between base stations adjacent to a terminal 
which is an object of the measurement, and the paths between the base stations 
5 adjacent to the terminal include non-line-of-sight multipaths. 

In step 350, the MULC 33 forms a path vector enumerating the adjusted 
time differences of the base station signals included into this path and determines 
its metric for each path of each generated set. 

In step 355, the MULC 33 selects a group of path vectors for each pair of 

10 base stations from the set of all possible path vectors of this pair of base stations, 
wherein the selected group of the path vectors contains each of the obtained 
adjusted time differences; wherein number of applications of each obtained 
adjusted time difference of the selected group of path vectors must not exceed 
number of applications of this adjusted time difference of any other group of path 

15 vectors obtained from the set of all possible path vectors, wherein values of path 
vector metrics of the selected group must not exceed values of path vector 
metrics of any other path vector group obtained from the set of all possible path 
vectors. 

In step 360, the MULC 33 forms weights of the adjusted time differences of 
20 signals of base stations for each pair of base stations using the selected group of 
path vectors of these base stations and the obtained accuracies of the adjusted 
time differences of signals of the base stations. For said all paths, weights are 
provided according to errors of the measured mutual time differences. 

In step 365, the MULC 33 determines a mutual time difference of signals 
25 of each pair of base stations as a weighted sum of all adjusted time differences of 
signals of base stations, wherein the weights of the adjusted time difference of 
signals of the base stations formed for the pair of base stations are used as 
weights. 

The claimed method of mutual time difference determination of base 
30 station signals in a cellular communication system has a number advantages 
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compared to prior inventions known to those skilled in the art. First, the claimed 
method enables the determination of the mutual time difference of signals of any 
pair of base stations regardless of the availability of the direct measurement of 
the adjusted time difference of signals of this pair of base stations. Secondly, the 
5 claimed method enables improved accuracy in determining the mutual time 
differences of base station signals. These advantages are achieved due to joint 
statistical processing of all adjusted time differences of base station signals. 

While the invention has been shown and described with reference to 
certain embodiments thereof, it will be understood by those skilled in the art that 
10 various changes in form and details may be made therein without departing from 
the spirit and scope of the invention as defined by the appended claims. 
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